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We have identified a class of zero-dimensional classical and quantum Heisenberg spin systems 
exhibiting anomalous behavior in an external magnetic field B similar to that found for the geo- 
metrically frustrated Kagome lattice of classical spins. Our calculations for the isotropic Heisenberg 
model show the emergence of a pronounced minimum in the differential susceptibility dM/dB at 
B sai t/3 as the temperature T is raised from K for structures based on corner-sharing triangles, 
specifically the octahedron, cuboctahedron, and icosidodecahedron. We also provide the first exper- 
imental evidence for this behavior dM/dB: It is exhibited by the giant Keplerate magnetic molecule 
{Mo72Fe3o} (Fe 3+ ions with spin s — 5/2 on the 30 vertices of an icosidodecahedron). The minimum 
in dM/dB is due to the fact that for low temperatures when B w B sa t/3 there exist two competing 
families of spin configurations of which one behaves magnetically "stiff" leading to a reduction of 
the differential susceptibility. 
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The magnetism of frustrated one-, two-, and three- 
dimensional lattice spin systems is a fascinating sub- 
ject due to the richness of phenomena that are ob- 
served [3, 0, 0- I n this Letter we report that one ef- 
fect of geometrical frustration, which so far has been 
reported [4| only for the theoretical model of classical 
spins on a Kagome lattice, already appears for a class 
of zero- dimensional materials, namely certain magnetic 
molecules hosting highly symmetric arrays of classical 
or quantum spins. These molecular units 5] contain a 
set of paramagnetic ions whose mutual interactions are 
described by isotropic Heisenberg exchange and where 
the mtermolecular magnetic interactions (dipole-dipole 
for the most part) are negligible as compared to in- 
tramolecular Heisenberg exchange. Magnetic molecules 
as zero-dimensional spin systems provide a new avenue 
for detailed exploration of the basic issues of geometric 
frustration. They are particularly appealing since they 
offer the prospect of being modeled unencumbered by 
some of the complications of bulk magnetic materials. 

We here report experimental and theoretical results 
for the occurrence of a striking anomaly in the differen- 
tial susceptibility dM/dB versus magnetic field B that 
is exhibited by the giant Keplerate magnetic molecule 
{Mo7 2 Fe 3 o} 00 This molecule features 30 Fe 3+ ions 
on the vertices of an icosidodecahedron that interact via 



nearest-neighbor (nn) isotropic antiferromagnetic (AF) 
exchange (J/ks = 1.57 K). Due to their near-perfect 
O^-symmetric coordination environment, the Fe 3+ ions 
represent ideal 5 = 5/2 spin centers with virually no 
single-ion anisotropy. We also present theoretical re- 
sults for the classical and quantum Heisenberg model 
showing that the same anomaly in dM/dB occurs for 
a class of geometrically frustrated zero-dimensional sys- 
tems, where spins mounted on the vertices of a triangle, 
octahedron, cuboctahedron, or an icosidodecahedron in- 
teract via nn isotropic AF exchange. As the temperature 
T is raised from K a deep narrow minimum in dM/dB 
emerges in the vicinity of one-third the saturation field 
B sat , which upon increasing T extends over a larger field 
interval and its sharp features progressively deteriorate. 
We attribute this phenomenon to a common topological 
property of these polytopes, namely that each is assem- 
bled from corner-sharing triangles. In the classical case 
the drop in dM/dB can be understood as a result of 
the interplay of two effects: In the immediate vicinity of 
^sat/3 a family of "up- up-down" (uud) spin configura- 
tions are energetically competitive with the continuous 
family of spin configurations of lowest energy 8] . How- 
ever, the uudspm configurations are magnetically "stiff", 
i.e. dM/dB w for low temperatures, and thus reduce 
the susceptibility of the system. 
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We write the AF Heisenberg Hamiltonian as 

ff = J^S m -S„+»B.^S„, (1) 

(m,n) n 

where J is a positive energy, the spin operators S n are in 
units of H, B is the external field, g is the spectroscopic 
splitting factor, \ib is the Bohr magneton, and (m, n) di- 
rects that the sum is over distinct nearest-neighbor pairs. 
The classical counterpart of Eq.Q is obtained by replac- 
ing each spin operator SL by y s(s + 1)5^ , where 5 n is 
a c- number unit vector 0, LU| • 

One very attractive feature of the polytopes under con- 
sideration is that their exact classical ground state energy 
is known For B < B sat it is given by 

E (B) = -^N A J C [1 + 3(JL^ ], (2 ) 

where J c = s(s + 1) J is called the classical Heisenberg 
exchange constant, B sat = 6J c //i c , /i c = g/i B y / s(s + 1), 
and TVa is the number of corner-sharing triangles (= 
4, 8, 20 for the octahedron, cuboctahedron, and icosido- 
decahedron, respectively). A plot of this quantity versus 
B/B sat is shown in Fig. [I] (solid curve). The ground state 
magnetic moment and differential susceptibilty are given 
by M (B) = -dE /dB and dM (B)/dB, respectively. 
For B = each spin system is decomposable into 3 sub- 
lattices of N/3 spins each; all spins of a given sublattice 
are mutually parallel; the sublattices are characterized 
by three coplanar unit vectors with angular spacings of 
120°. The magnetization of the system is linear in B un- 
til .Bgat and constant (fully saturated configuration) for 
larger fields. The linear rise with B can be pictured in 
terms of the folding of an "umbrella" 8] defined by the 
three sublattice unit vectors as they close towards the 
field vector B. 
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FIG. 1: Total energy vs. magnetic field for T — K for the 
classical AF triangle, octahedron, cuboctahedron, and icosi- 
dodecahedron. The solid curve is given by Eq.|5| The dashed 
curves correspond to collinear structures discussed in the text. 



Also shown in Fig. are the energy curves for three 
other specific configurations of interest. These are con- 
figurations where the three (unit) spin vectors associated 
with each triangle are constrained to be collinear and 
the resultant vector is either parallel or anti-parallel to 
B. These configurations are labeled as uuu (up-up-up), 
uud, and ddu. For each of these collinear configurations 
the magnetic moment of the polytope is independent 
of B and thus dM/dB vanishes and one can describe 
these configurations as being magnetically "stiff". The 
fully saturated uuu configuration is of minimal energy 
for B > B sat . The uud configuration is of special inter- 
est since its energy coincides with the minimal energy of 
the spin system for B = B sat /3 and exceeds the minimal 
energy for any other choice of field. For T — K and 
for any choice of B other than B sat the uud configura- 
tion will not play a role. However, for T > K and for 
B in the vicinity of B sat /3 a significant contribution to 
the partition function will arise from the set of configu- 
rations derived by infinitesimal modifications of the uud 
configuration. These slightly modified uud configurations 
lead to a reduction of the differential susceptibility of the 
system because of their magnetic stiffness. Our quali- 
tative considerations for T > K are confirmed by the 
results of our classical Monte Carlo simulations for the 
three polytopes as shown in Fig. [2j Fig. 03 displays the 
results for a classical model of {Mo72Fe3o}, namely 30 
classical spins on the vertices of an icosidodecahedron, 
as substitutes for quantum spins with 8 = 5/2. As T 
is increased from K a sharp narrow drop emerges that 
is situated at B sat /3 (vertical dashed line). As T con- 
tinues to increase the drop extends over a larger interval 
and its sharp features progressively wash away. One also 
observes a temperature dependence of the field associ- 
ated with the minimum in dM/dB, i.e., it decreases with 
increasing T. 
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FIG. 2: Low-temperature (kBT/J c = 10 -2 ) simulational re- 
sults for dM/dB vs. B for classical spins on the octahedron, 
cuboctahedron, and icosidodecahedron. 

The relevance of these theoretical results to real mag- 
netic materials is demonstrated by our experimental data 
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for the differential susceptibilty of the giant Keplerate 
magnetic molecule {Mo72Fe3o}. The magnetization was 
measured at 0.42 K in pulsed magnetic fields up to 23 
Tesla (sweep rate 15000 Tesla/s) at the Okayama High 
Magnetic Field Laboratory by using a standard induc- 
tive method. The sample was immersed in liquid 3 He to 
maintain good contact with the thermal bath. The exper- 
imental results for dM/dB (in arbitrary units) are shown 
in Fig. 01 and the drop at about B sat /3 is clearly evident. 
However, the data resembles the simulational curve for 
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FIG. 3: Differential susceptibility dM/dB versus B for the 
classical Heisenberg model of {Mo72Fe3o} obtained by Monte 
Carlo simulations for temperatures given in the legend. 



2 K (see inset Fig. 0}, not 0.42 K, perhaps suggesting 
an elevated effective spin temperature due to the high 
sweep rate. To clarify this point we also measured the 
magnetization in steady fields by a capacitance method 
in the range up to 7 Tesla for T = 16 mK and 0.73 K 
and obtained close agreement with the pulsed-field data. 
The considerable broadening of the drop in dM/dB may 
point to the occurence of a staggered field jl2| , since the 
principal axes of the Fe ions are not strictly equivalent. 
A specific suggestion [T^j is that the broadening is due to 
Dzyaloshinskii-Moriya terms supplementing the isotropic 
Heisenberg model, originating from possible low symme- 
try of the nearest-neighbor Fe-Fe bond. Further study of 
this issue is warranted. 

To explore the role of quantum effects we have calcu- 
lated dM/dB for the triangle and the octahedron of spins 
with arbitrary s as well as for a cuboctahedron (TV = 12) 
with s = 1/2 and s = 1. For the latter system this in- 
volves numerical diagonalization of matrices defined on 
a Hilbert space of dimension 3 12 (= 531441). Even by 
fully exploiting the symmetries of the Hamiltonian this 
is at the limit of present day computing capabilities. The 
results for the cuboctahedron with s — 1 are shown in 
Fig. for different temperatures. As in the previous fig- 
ures one again encounters a strong reduction of dM/dB. 
The minimum is located exactly at M sat /3, but since B 
and M are not strictly proportional for a quantum sys- 
tem the drop occurs for fields slightly larger than B sat /3. 
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FIG. 4: Experimental results (in arbitrary units) for a sample 
of {Mo72Fe3o} performed at 0.42 K using a pulsed-field tech- 
nique. In the inset Monte Carlo results for 0.42 K and 2.0 K 
are given. 



For T = K M vs. B can be described as a "stair- 
case" of 12 steps originating from ground state Zeeman 
level crossings and dM/dB consists of a set of Dirac delta 
functions at the crossing fields. For T > K the abrupt 
magnetization steps are smoothed out and dM/dB ex- 
hibits finite peaks. Our results for the triangle and the 
octahedron for general spins s exhibit the same overall 
behavior seen in Fig. 

One can understand that the pronounced minimum in 
the susceptibility occurs for classical as well quantum 
spins by examining the partition function for the par- 
ticularly simple example of the triangle where the results 
can be obtained by exact analytical methods. For integer 
spins s the quantum partition function may be written 
as 

3s 2 

Z(t, b) = (sinh(6cro)) _1 ^ G n sinh(&<j n ), (3) 
wher e b = f i c B/(k B T), t = k B T/J ( 

\)/y/ 8 {8+l), G n = r n /0>(*+l) 

tiplicity factor, namely the number of distinct ways of 
achieving total spin n upon adding three distinct (inte- 
ger) quantum spins s. In particular T n = 2(n + -|) for 
< n < s and T n = 3(s + \) - (n+ \) for 5 + 1 < n < 3s. 
The analogous formulas are easily derived for half-integer 
spins s. Formula (0) for Z is very similar to that for the 
classical Heisenberg triangle which may be written as [Icj 



cr n = (n + 
and T n is the mul- 



Z(t,b) 



JO 



dSG(S)e-^smh(bS). (4) 



Here G(S) = 2S for < S < 1 and G{S) = 3 - S for 
1 < S < 3, arising from considering the geometrical vol- 
ume available to three unit vectors such that the magni- 
tude of their vector sum lies within a shell of radius S and 
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unit thickness. Indeed it is straightforward to verify that 
in the limit s — > oo the quantum result (Eq. JJJJ)) agrees 
with the classical formula (Eq. Q). In the quantum for- 
mula the multiplicity factor corresponds to the classical 
geometrical function G(S). Each of these quantities has 
two distinct branches, depending on whether n is in the 
interval [0, s] or [s + 1,3s] or whether S is in the inter- 
val [0,1] or [1,3]. In fact, the existence of two distinct 
branches becomes manifest in various higher derivatives 
of Z(t, b) at nonzero temperatures for fields in the vicin- 
ity of B = B sat /3. For < t <C 1 there exists a narrow 
field range at about B sat / 3 such that each of the functions 
exp(-a 2 L /(2t))smh(ba n ) and exp(- S 2 / (2t)) sinh(6S) has 
a very narrow maximum for <j n « 1 and S « 1 but nev- 
ertheless samples the two branches. This is the mathe- 
matical orgin of the pronounced minimum in dM/dB at 
B = B sat /3. 




FIG. 5: Differential susceptibility dM/dB versus B / B sa t of 
the quantum Heisenberg cuboctahedron (s — 1) for values of 
UbT I J shown in the legend. 



Plateau like structures in the magnetization versus B 
in various two- and three-dimensional lattices built of 
corner-sharing triangles lattices at one-third of the sat- 
urated moment have been under investigation for the 
past two decades as an expression of geometric frustra- 
tion [J, 0, 0, 0, E3 Moreover, theoretical studies of 
the classical Heisenberg antiferromagnet on the Kagome 
lattice show that dM/dB has a pronounced minimum at 
one-third of B sat |4(. However, the study of selective mag- 
netic molecules such as {Mo72Fe3o} can give new insights 
for this subject since such molecules are much better ac- 
cessible both experimentally and theoretically. 

In summary, we have shown that for a class of ge- 
ometrically frustrated magnetic polytopes, namely the 
octahedron, the cuboctahedron and the icosidodecahe- 
dron, field-induced competitive spin configurations exist 
which manifest themselves in an pronounced minimum 
in the differential susceptibility dM/dB in the vicinity of 
Bs&t/3- We have also reported the first experimental ob- 
servation of this effect. Our data for the giant Keplerate 



magnetic molecule {Mo72Fe3o} are consistent with our 
classical Monte Carlo results for the icosidodecahedron. 
Furthermore, we have shown that this feature reflects 
a general intrinsic property of the very building block 
of these specific polytopes, namely the simple AF equi- 
lateral Heisenberg spin triangle, and emerges for both 
classical and quantum spins. Moreover, our theoretical 
calculations for each of these polytopes show that the 
specific heat versus B also exhibits anomalous behavior 
in the vicinity of B sat /3 [lfij |. A measurement of this 
quantity for {Mo72Fe3o} at very low temperatures would 
be of great interest. 
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